Dynamics of conduction blocks in a model of paced cardiac tissue 
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We study numerically the dynamics of conduction blocks using a detailed electrophysiological 
model. We find that this dynamics depends critically on the size of the paced region. Small pacing 
regions lead to stationary conduction blocks while larger pacing regions can lead to conduction blocks 
that travel periodically towards the pacing region. We show that this size-dependence dynamics can 
lead to a novel arrhythmogenic mechanism. Furthermore, we show that the essential phenomena 
can be captured in a much simpler coupled-map model. 
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I. INTRODUCTION 

The coupling of the electrical excitation to the contrac- 
tile forces in the heart is essential to the blood supply 
to the body. Abnormal electrical activity, in particular 
arrhythmias, can have dire consequences. The most se- 
rious of all, ventricular fibrillation (VF), is characterized 
by disordered electrical activity within the ventricles and 
leads to death within minutes. 

Unfortunately, the precise cause for the onset and 
maintenance of VF remains elusive. It is believed, how- 
ever, that electrical wave reentry plays an important role. 
Reentry is initiated when an electrical wave is locally 
blocked, leading to a broken wave front that can re-excite 
the tissue behind the conduction block. There are several 
ways to create a conduction block. Perhaps the most in- 
tuitive one is an anatomical tissue heterogeneity through 
which the electrical wave fails to propagatepj Q- Al- 
though these heterogeneities exist and play a role in the 
initiation of VF, recently a new mechanism for conduc- 
tion block in perfectly homogeneous tissue has been de- 
scribed in models [1,0,01 and experiments [a, El- Key 
contributor to this dynamical heterogeneity is electrical 
alternans which is characterized by a beat to beat os- 
cillation in the action potential duration (APD, defined 
below) under rapid pacing conditions. Clinically, the oc- 
currence of T wave alternans in ECGs has been associ- 
ated with sudden cardiac arrest [1, @ • In isolated cells, 
the onset of alternans can be determined from the resti- 
tution curve which relates the APD to the diastolic inter- 
val (DI), the time interval between two successive action 
potential. If this curve has a slope larger than one, it 
is easy to see that a period doubling bifurcation devel- 
ops, resulting in an APD that is alternatingly short and 
long [l(3,[ll|. This simple vision, however, is most likely 
not complete as memory effects |l2j and calcium cycling 
effects [3, 0| have been shown to play a role in the 
mechanism leading to alternans. 



In the case of spatially extended systems, the possi- 
ble dynamical behavior becomes more complicated. The 
conduction velocity also depends on the DI, which can 
lead to a spatial modulation of the alternans, called dis- 
cordant alternans. 0, Q. During discordant alternans, 
the APD is following a long-short-long pattern in one 
region of the cable and a short-long-short pattern in an- 
other. Separating these regions are nodes where the APD 
is constant. They can be either stationary or traveling. 
In addition, during discordant alternans, the amplitude 
of the alternans (the difference between APDs in sub- 
sequent beats) can grow when moving away from the 
pacing site. Hence, at a critical distance, provided the 
amplitude of alternans is large enough, the tissue can no 
longer support a traveling wave for each excitation and a 
diffusive conduction block[l5l| develops: a wave will prop- 
agate over a given distance along the cable and will then 
come to a stop. 

Here, we study the dynamics of the conduction blocks 
using a detailed electrophysiological model and a pac- 
ing domain of variable size. We present a quantitative 
and qualitative description of the observed patterns and 
present results from numerical experiments in cables and 
sheets of cardiac tissue. Our main results are: (i) in ad- 
dition to a stationary conduction block, rapid pacing can 
lead to periodic patterns of moving conduction blocks, 
(ii) the type of conduction block is critically dependent 
on the size of the pacing region, and (iii) these effects 
combined can lead to a novel mechanism for reentry. 
Furthermore, we show that a simple coupled-map model 
gives qualitatively similar results. 



II. MODELING OF A CABLE OF CARDIAC 
CELLS AND TERMINOLOGY 

To describe the electrophysiological properties of car- 
diac tissue we used the reaction-diffusion equation: 



dtV = DS7 2 V - (£ wn I lon + I sUm )/C 



(1) 
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where V is the transmembrane potential, C = 1 /zF cm 2 
is the membrane capacitance, D V 2 , with D = 



2 



0.001cm 2 s _1 , expresses the inter-cellular coupling, Ii on 
represents the different trans-membrane currents and 
Istim is the applied pacing current. The ionic currents 
in Ii on are governed by nonlinear evolution equations 
coupled to V and here we have chosen to use a modi- 
fied version of the Luo-Rudy model ^(| to describe these 
currents. A detailed description of this model and our 
modifications are given in appendix^ The integration 
scheme we used was forward Euler with space and time 
discretizations of Sx — 0.025 cm and 5t = 0.02 ms, re- 
spectively. Finally, the space constant was determined to 
be 0.06 cm. 

The cable was paced at one end with a period T 
through a 1 ms long constant current stimulus of I s um — 
— 80 fi A/ cm 2 |2g The pacing domain consisted of n grid- 
points, where n was varied between 1 and 20, thus corre- 
sponding to a pacing domain ranging in size from 0.025 
cm to 0.5 cm 27] As we will see below, changing the 
number of stimulated grid points can have a dramatic 
effect. On the other hand, changing the size of the simu- 
lation cable from 10 cm, the length used here, to 5 cm did 
not reveal any qualitative differences. This, along with 
the fact that the propagation speed from endocardium 
to epicardium can be much smaller (about 17 cm/slpT) 
than the one used here (about 51 cm/s), implies that the 
observations reported here might be relevant to human 
ventricles. Finally, additional simulations performed us- 
ing either a stronger stimulus or a longer stimulus showed 
little effect. 

Wc define the beginning of an action potential at po- 
sition x as the moment when V(x,t) crosses a threshold 
value Vth from below. Similarly, the end of the action po- 
tential is defined as the moment when V(x,t) decreases 
below Vth- The APD is then defined as the time differ- 
ence between these two events (see Fig. 1). Then, the 
DI is defined as the time between the end of one action 
potential and the beginning of the next (see Fig. 1). A 
conduction block for the m th stimulus occurs at the i th 
grid point when an action potential was elicited at the 
jth _ ^ p i n t by the m th stimulus but not at the i th 
one. Vth was chosen to be equal to -60 mV and changes 
in the value of Vth (between -80 mV and -50 mV) did 
not affect results significantly. In addition, it is useful to 
define a waveback and a wavefront. A wavefront is the 
boundary between a region at rest and an excited region 
when the latter region is invading the former while the 
waveback is the boundary when the former is invading 
the latter. 

The pacing periods used here (« 160 ms) are short 
compared to the normal pacing period in humans. Of 
course, this pacing period is dictated by the the choice of 
the electrophysiological model used here. 



III. RESULTS IN A CABLE 

Let us first consider a one dimensional cable. For small 
enough pacing period, a period doubling instability to al- 
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FIG. 1: top: Schematic time course of the transmembrane po- 
tential at a point during normal cardiac rhythm (i.e. without 
alternans). The diastolic interval (DI) is the period of time 
during which the transmembrane potential is below a given 
value and the action potential duration (APD) is the period of 
time during which the action potential is over this threshold 
value. Bottom: Schematic drawing of the APD along the ca- 
ble during an even beat (solid line) and an odd beat (dashed 
line) in the discordant alternans regime. In the case of normal 
cardiac rhythm (i.e. without alternans), both lines would be 
the same and the APD would be constant in space (except 
at the pacing site (resp. the end of the cable) where bound- 
ary effects tend to lengthen (resp. shorten) the APD). In 
the case of concordant alternans, the APD during successive 
beats would be different but constant along the cable (apart 
from boundary effects). 



ternans takes place. After a transient regime where the 
amplitude of alternans grows, a steady state is reached. 
This regime is characterized by the fact that each stimu- 
lus will elicit an action potential all along the cable and 
that the action potentials elicited by two successive stim- 
ulus at a given point in space will have different durations 
(a short action potential follows a long one which follows 
a short one and so on)[l8j. This allows us to define, 
at a given point in space x, the amplitude of alternans, 
a(x, m), for the m th stimulus as the difference between 
the APD due to this stimulus and the APD due to the 
previous stimulus multiplied by (—1)"' so that a does 
not change sign every stimulus. In the case of concor- 
dant alternans, a is roughly constant in space and does 
not change sign. In our simulations, however, since the 
size of the simulated cable is big ger than the character- 
istic wavelength of alternans [Ljfl , we observe discordant 
alternans. In this case a(x, m) is not constant in space 
and changes sign when going along the cable. For ex- 
ample, if a point in the cable with positive a exhibits a 
sequence of beats that is long short long short... then a 
point of the cable with negative a will exhibit a sequence 
of beats that is short long short long.... Consequently, be- 
tween regions of opposite a, there is a point where a = 0, 
corresponding to an alternans node. 

In the case of discordant alternans, one can observe two 
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FIG. 2: al and a2: Grey scale time-space plot of the alternans 
amplitude |o(a;,m)| with black corresponding to a(x,m) = 
and white corresponding to large values of the amplitude. 
The pacing period in both figures is 162 ms and the number 
of pacing grid points is n = 6 in al and n — 10 in a2. In 
al, one can see discordant alternans where the nodes, regions 
of zero amplitude (and thus in black), are not moving. In 
a2, the discordant alternans leads to nodes that are traveling 
towards the pacing site, bl to b5: successive profiles of V(x), 
taken at 20 ms intervals, during a conduction block event. 
The waves are initiated on the left hand side of the cable. 
In bl, one can see an action potential propagating from the 
left towards the right with its typical sharp wavefront. The 
smoother waveback of the wave due to the previous stimulus 
is on the right hand side of the plot. The wavefront propa- 
gates faster than the wave back (bl to b3) and finally reaches 
a region that has not yet recovered and can not be rapidly ac- 
tivated. Reaching this region (b4), the wave comes to a stop 
and fails to propagate further (this is called throughout this 
paper a conduction-block). In b5, one can see that the sharp 
wavefront has disappeared. The cable is 10 cm long and V 
ranges from -100 mV to 100 mV. 



distinct regimes (see [13 for a full theoretical analysis) 
that are characterized by the spatio-temporal evolution 
of a. The first regime, called standing alternans, is char- 
acterized by stationary alternans nodes. This is shown 
shown in Fig. |2]al, where we show a gray-scale space- 
time plot of a. In contrast, the alternans nodes in the 
second regime, called traveling wave alternans, are non- 
stationary and are traveling towards the pacing site. This 
is illustrated in Fig|3a2 which shows that a is changing 
slowly in time with a given period. In addition, in this 
regime, a = at the pacing site and \a\ increases when 
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FIG. 3: Position of the conduction block as a function of 
the stimulus number m for different pacing period T (in ms) 
and pacing region sizes n. In a, the inset shows that the 
conduction block leads to a 2:1 rhythm. In both b and d, 
the conduction block travels towards the pacing site during a 
number of stimuli, followed by a small number of stimuli dur- 
ing which there is no conduction block. In b and d the thick 
lines correspond to the stimuli where no conduction block is 
observed. One should note the similarity between the trajec- 
tories of conduction blocks for T = 159 ms and n = 10 (resp. 
n = 6) and the trajectories of the alternans nodes for T — 162 
ms and n = 10 (resp. n — 6) shown in figurc|5|al and a2. 



moving away from the pacing site. 

A change in the dynamics of the alternans nodes was 
observed earlier in the theoretical study of Echebarria 
and Karma |l9j . There, this change was accomplished 
by changing the properties of the tissue (i.e. by changing 
the model). Here, in contrast, the only difference between 
the simulations shown in Fig|2|al and a2 is the size of the 
pacing region: n = 6 in al and n = 10 in a2. This finding 
already shows that the size of the pacing region may play 
a critical role in the alternans induced arrhythmias. A 
rationale for the effect of changing the size of the pacing 
region will be given later in the manuscript (end of Sec. 
V). 

Decreasing the stimulation period further leads to con- 
duction blocks. In other words, some stimuli are not able 
to create an action potential that propagates all along 
the cable. Two different types of conduction blocks were 
observed: 1) conduction block at the pacing site; the stim- 
ulus is not able to elicit an action potential at all and 2) 
conduction block away from the pacing site; the stimu- 
lus creates an action potential that begins to propagate 
along the cable and this action potential fails to propa- 
gate at a given point in space away from the pacing site. 
The latter type of conduction block is illustrated in Fig. 
[21 bl-5. Hence, a conduction block can be characterized 
by the point in space where it took place and also by the 
index of the stimulus that created the action potential 
which failed to propagate at this point. 
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In our simulations, we found that the dynamics of the 
conduction blocks depends on both the size of the paced 
region and on the pacing period. Examples of the ob- 
served dynamics of the conduction blocks are illustrated 
in Fig. |31 where we plot the position of the conduction 
block as a function of the pacing cycle number m. In 
fig. Oil and c, the block occurs at a fixed location (ei- 
ther away from the pacing site or at the pacing site) 
once every other stimulus (i.e. a 2:1 rhythm). Exam- 
ples of more complex behavior are illustrated in Fig. [jJj 
and d where the location of the conduction block is non- 
stationary and forms a periodic pattern. In these cases, 
a conduction block forms at a certain location away from 
the pacing site. In subsequent beats, the position of this 
conduction block moves closer to the pacing site (during 
this regime every other stimulus is blocked). Eventually, 
the block disappears, after which no conduction blocks 
are observed during several pacing cycles (i.e. propa- 
gation with a discordant alternans pattern). Then, the 
conduction block reforms at precisely the same location 
as previously and the sequence restarts. 

This gives rise to a periodic pattern of conduction 
blocks with a certain "amplitude" , defined here as the 
distance between the site where the conduction block first 
forms and the site where the conduction blocks disap- 
pears. We found that this amplitude can take on values 
between 2 cm (see Fig. |3|i) and 5 cm (see Fig. Eb) and 
that for a given pacing period, both the amplitude of the 
conduction block wave and its period were dependent on 
the number of pacing grid points. We also found that the 
trajectory of the conduction blocks in Fig. EJb), which is 
close to the onset of conduction blocks, is strikingly simi- 
lar to the dynamic pattern of the location of the alternans 
nodes shown in Fig. [5] a2. These alternans nodes also 
display periodic movement towards the pacing site with 
slow movement followed by fast movement. This move- 
ment is presumably due to the use of zero flux boundary 
conditions which favors the existence of an extremum at 
the boundaries: slow movement corresponds to an ex- 
tremum of the alternans amplitude near the boundary 
while fast movement corresponds to an alternans node 
near the boundary. A phase diagram representing the 
nature of conduction blocks as a function of the number 
of pacing grid points and of the period is shown in Fig. 

It is also worth mentioning that the use of slightly ir- 
regular pacing times, modeled via the inclusion of a noise 
term with a variance of 1 ms, did not change significantly 
the phase diagram. 



IV. RESULTS IN A SHEET OF CARDIAC 
TISSUE 

From the above, it should now be clear how a slight 
inhomogeneity in the size of the pacing area can lead to 
reentry. Envision a 2D sheet of cardiac tissue (modeled 
by Eq^ using two space variables instead of one), paced 
from one side by a domain containing two widths, n for 



the upper part and n* for the lower part, as shown in Fig. 
13 Neglecting spatial coupling in the direction perpendic- 
ular of the wave propagation, the lower and upper part of 
the tissue will see waves originating from domains with 
different sizes. Reentry should be possible when the spa- 
tial locations of the conduction blocks in the two domains 
are significantly different. Then, the electrical wave will 
be blocked in one part of the tissue while propagating 
normally in the other part, leading to the excitation of 
the tissue behind the conduction block and to reentry. 



We can estimate the likelihood of reentry in 2D, paced 
by domains of size n and n* , based on our ID results. 
Reentry is likely when (i) a cable paced with n shows a 
conduction block while a cable paced with n* does not 
exhibit a conduction block, or vice-versa, or (ii) the dif- 
ference in spatial location between the conduction blocks 
in cables paced with n and n* becomes large. We found 
that for n* = nil, one or both of these conditions 
were met in most of the phase space for which conduc- 
tion blocks occur away from the pacing site (i.e. filled 
circles in Fig. Q). In most cases, the conditions were met 
within a few hundred stimuli while occasionally a larger 
number of stimuli was necessary. 



Of course, the preceding arguments neglect spatial cou- 
pling perpendicular to the wave propagation direction. 
Nevertheless, we found that the phase diagram was a 
valuable predictor for the occurrence of reentry in full 
scale 2D simulations. A typical example of our 2D sim- 
ulations, performed using the same numerical scheme as 
in our one dimensional simulations, is displayed in Fig. 
[5j where we show a series of gray-scale plots of the mem- 
brane potential, with black corresponding to repolarized 
tissue and white corresponding to depolarized tissue. In 
this example, the width of the two pacing domains dif- 
fered by a single grid point: n* = 8 vs. n = 9. Both 
domains are paced with a constant period of 158 ms, 
leading to traveling conduction blocks with slightly dif- 
ferent periodicity. During the roughly 100 first stimuli, 
the position of the conduction block was nearly identical 
in the upper and lower parts of the tissue. During the 
subsequent stimulus, however, a wave block remains in 
the upper part of the tissue but not in the lower part 
of the tissue, which now allows unblocked propagation 
(see Fig. |3Jd and d). As time progresses, the propagating 
wave in the lower part reenters the upper part behind 
the wave block, leading to a spiral wave. For clarity, we 
have stopped stimulating the tissue once the reentry ap- 
peared, although we have verified that continuous stim- 
ulation also resulted in long-lasting reentry. In addition, 
we have checked that smaller domains were also able to 
produce reentry and that using a line of stimulation with 
a thickness varying randomly can also lead to sustained 
reentry. 
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FIG. 4: Phase diagram showing the type of conduction 
block: filled circles represent periodic conduction block waves, 
crosses represent stationary conduction blocks at the pacing 
site, triangles represent stationary conduction blocks away 
from the pacing site and filled boxes represent no conduc- 
tion block (i.e. discordant or concordant alternans.). T is 
expressed in ms. 



V. COUPLED MAP MODEL 

The results from the full ionic model can be repro- 
duced by a coupled map model similar to the ones used 
in previous work 0, Ei- To describe the dynamics of 
alternans nodes and conduction blocks, we consider the 
following two fields: T up (x, n), which is the time at which 
an action potential is elicited in x due to the n th stimulus 
and T rep (x,n), which is the time at which it ends. The 
equation for the repolarization time T rep is the same for 
both the cable and the pacing domain: 



T rep (x, n) = T up (x, n) + APD(DI(x, n)) + 

+ ( 2 V 2 T rep (x, n) + wd x T rep (x, n) (2) 

The first two terms in this equation express the fact that 
a single cell repolarizes at time T up plus the duration of 
the action potential. The latter is taken from the restitu- 
tion curve APD(DI). For the paced domain, this curve 
is calculated using a single cell while for the cable it is cal- 
culated 1 cm away from the pacing domain, taken to be 
large. In this region of the cable, we found that the resti- 
tution curve was minimally affected by the proximity of 
the pacing domain and that the dispersion of repolariza- 
tion due to alternans was minimal. The two restitution 
curves are shown in Fig. [SJ1. The last two terms take 
into account spatial effects: the diffusion term indicates 
the perturbation due to inter-cellular coupling (we con- 
sider here that the repolarization wave is a phase wave 
and does not correspond to the propagation of a wave 
back) and the last term expresses the asymmetry intro- 
duced by the pacing at one end of the cable [l9j. The 
length scales £ and w were chosen here to reproduce the 
phase-diagram of Fig. 0]in a semi-quantitative fashion. 
The equation describing T up in the pacing domain is 
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FIG. 5: Series of gray-scale plots showing the the birth of 
a spiral wave. Time (shown in ms) is arbitrarily set to at 
the start of the last stimulus. The gray-scale represents the 
membrane potential values with white corresponding to max- 
imum depolarization and black corresponding to maximum 
repolarization. At t=0 ms, one can see the wave elicited by 
the previous stimulus propagating from left to right. Because 
of the discordant alternans, the APD differs strongly across 
the sheet with the left most clear stripe corresponding to a 
region where the APD is long (it is not a propagating wave 
despite the cells in this region are depolarized |lal ) and the 
central dark stripe corresponding to a region where the APD 
is short. At t=100 ms, one can see on the left hand side of 
the plot the wave elicited by the next stimulus while on the 
right hand side one can see the waveback of the wave seen 
in the t=0 ms frame. At t=200 ms, a conduction block oc- 
curred in the top part of the sheet while in the bottom part 
the wave was able to continue to propagate. This partial wave 
block results in a broken wave front which leads to the birth 
of reentry as can be seen in the t—240 ms, 280 ms and 320 
ms frames. As this wave front curls and attempts to reenter 
previously excited tissue with a long repolarization time (gray 
stripe in t=240 ms frame), it breaks up into two spiral tips 
(t=280 ms). Finally at t=760 ms, one can see further insta- 
bilities have led to a disordered activity similar to ventricular 
fibrillation. The geometry of the computational domain is 
also shown schematically where the gray region, exaggerated 
for clarity, corresponds to the cells that are paced. The width 
of the bottom part of this region is n* = 8 and of the top part 
n = 9 grid points. 



in the cable it is given by 

T up (x, n) = T n + I dx' 



o c(DI(x',n)) 



(3) 



where c is the propagation speed of the wave front. This 
propagation speed is a function of the DI, which it- 
self is of course a function of x and is the difference 
between the arrival of the n th action potential at x 
and the repolarization following the (n — l) th stimulus: 



T up (x 1 n) = T„, with T n being the pacing period, while DI(x,n) = T up (x,n) — T rep (x,n — 1). 



The restitution curve also determined the occurrence of 
a conduction block: we defined a conduction block to take 
place when the restitution curve had no value for the dias- 
tolic interval computed from Eq. [3J In this case, we com- 
puted the repolarization times between the origin and the 
position of the conduction block using Eq. along with 
boundary condition T rep = T up at the site of the conduc- 
tion block. The repolarization time for the remainder of 
the cable (i.e. behind the conduction block) was set to 
the repolarization time computed at the previous stimu- 
lus, except in a small transition region, in size equal to £, 
immediately behind the conduction block. There, a sim- 
ple sigmoidal function was used that smoothly connected 
the repolarization times on both sides of this region. 

The reduced model is able to capture the essential fea- 
tures of the full ionic model, including the striking de- 
pendence on the size of the pacing region for both the al- 
ternans nodes and the conduction blocks dynamics. The 
latter is shown in Fig. [^,-c where we plot the position 
of the conduction block as a function of pacing cycle for 
w = 0.025 and £ 2 = 0.04 which are comparable to the 
values presented in Table I of Ref. [Tjj when consider- 
ing a full ionic model. All qualitative features of Fig. |3| 
can be readily recognized, including traveling conduction 
blocks with a large amplitude (a) j28j, traveling conduc- 
tion blocks with a small amplitude (c), and stationary 
conduction block away from the pacing domain (b) . Not 
shown here, but also found within the reduced model, is 
2:1 conduction block at the pacing site. For the values 
of w and £ employed here, the boundaries between the 
different types of conduction block in the phase diagram 
derived from the full model and the coupled-map model 
differed by at most 5 ms. 

Both the full ionic model and coupled map model show, 
for certain pacing periods, a transition from stationary 
conduction blocks to traveling conduction blocks as the 
size of the pacing domain is increased. Interestingly, a 
similar transition also occurs in the dynamics of alter- 
nans nodes when one increases this size. To understand 
the effect of changing the size of the pacing region on 
the nature of alternans we consider the equation for the 
amplitude a of the alternans derived in 



cqu; 
.19] = 



dta = aa — \ dx'— — wd x a + £, 2 d xx a — ga 3 , (4) 



A 



Here, aa represents a linear growth term, ga 3 is a non- 
linear restabilizing term and A, w and £ are length scales 
which are a priori different from the ones w and £ used in 
Eqs. [3 El The transition between traveling nodes in our 
simulations (large pacing domains) and standing nodes 
(smaller pacing domains) can be understood when real- 
izing the size of the pacing domain dictates the boundary 
condition at x = of Eq. Pacing a cable using a large 
piece of tissue results in a pacing domain that essentially 
behaves as a single isolated cell. If the restitution curve 
of a single cell (solid line in Fig. Eh>) has a slope less than 
one, this domain will not display alternans and the am- 
plitude equation needs to be solved with boundary condi- 
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FIG. 6: a-c: Position of the conduction block as a function 
of the stimulus for different pacing period and pacing region 
sizes calculated using the coupled-map model (a: T=161 ms, 
n=20, b: T=161 ms, n=10, c: T=156 ms, n=12). In a and c 
the thick lines correspond to the stimuli where no conduction 
block is observed, d: The restitution curve used in the model 
for the single element (solid line) and for the cable (dashed 
line). Note that the slope of the single cell curve is smaller 
than the slope of the whole cable curve. 



tion a = 0. For a smaller pacing domain, however, spatial 
coupling becomes important and the relevant restitution 
curve becomes steeper (dashed line in Fig. Eh>), which 
can lead to alternans. Thus, for smaller domains the 
relevant boundary condition for the amplitude equation 
is d x a = 0. A standing wave solution to the amplitude 
equation can be described by a — aocos{2 j n Xx + (f) , where 
1/A is the separation between the nodes. This solution 
(with do ^ 0) can only satisfy the d x a = boundary 
condition and not the a = boundary condition. This 
is easy to see when one considers the standing wave so- 
lution in Eq. 0] at x = 0: the only non zero term at 
the r.h.s. is d x a which requires oqItt/X = and finally 
a.Q = 0. Hence, a standing wave solution is only possible 
for small pacing sizes and increasing the size of the pac- 
ing domain suppresses the standing wave solution. It is 
likely that this qualitative change in the nature of allowed 
solutions underlies the transition between stationary to 
traveling blocks observed in the ionic and coupled map 
model. 



VI. DISCUSSION 

To conclude we have found that the spatio-temporal 
structure of alternans and the dynamics of conduction 
blocks is strongly influenced by both the size of the paced 
domain and the pacing period. This dependence provides 
a novel arrhythmogenic mechanism which we illustrate in 
a homogeneous two dimensional sheet of tissue, paced by 
two domains that vary slightly in size. The reentry is ini- 
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tiated through pacing with a constant period and hence 
does not require an abruptly changing pacing frequency 
nor a symmetry-breaking change of location as in previ- 
ous studies 0. We have found, using different ionic mod- 
els, including a detailed canine model H an d the simpli- 
fied three variable Fenton-Karma model |20Ll2l|. that the 
results we present here are partially model-dependent. 
This is perhaps not surprising, as other arrhythmogenic 
mechanisms have been show to depend on the details of 
the electrophysiological model |22|. Nevertheless, a re- 
duced coupled map model was able to reproduce, both 
qualitatively and, to some degree, quantitatively, the re- 
sults presented here. Consequently, we expect that the 
dynamics of conduction blocks will depend critically on 
the pacing protocol in a wide range of detailed models. 

Of course, the ultimate test should come in the form 
of experimental studies similar to the one conducted by 
Fox et al but with a varying in size stimulation region. 
It should be possible to conduct quantitative studies of 
Purkinje fibers, which conduct the electrical stimulus in 
an actual heart to the ventricles and which penetrate 
the heart wall to varying depths. These fibers can be 
isolated, resulting in linear strands of cardiac tissue 0. 
On the theoretical side, a further extension of this work 
will be the formulation of a coupled map model which 
takes into account the effects of coupling transversely to 
the propagation direction of the wave. 
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APPENDIX A: DESCRIPTION OF THE MODEL 



The Luo-Rudy 
ments (see e.g. ,2. 



model ^(|, and its subsequent refine- 
J), have been widely used either in their 



original form or in modified forms in numerical studies of 
wave propagation in cardiac tissue [24|- The model de- 
scribes the voltage and time dependence of the ionic cur- 
rents Tii on Ii on used in equation for the transmembrane 
potential V of a single cell: 



C 



dV 
~dt 



V ■ T 



(Al) 



where C is the membrane capacitance. In the original 
Luo-Rudy model, the total ionic current is given as 



-(iNa + Isi + Ik + Iki + 1 



Kp 



(A2) 



where I^ a — GN a m 3 hj(V — £^v a ) is the fast sodium 
current, I S i = G S idf(V — E S i) is the slow inward cur- 
rent representing the L-type calcium current, Ik = 
GkxX\ (V — Ek) is the time dependent potassium cur- 
rent, Iki = GkiK^oo{V — Eki) is the time indepen- 
dent potassium current, Ik p = Gk p K p (V — Ek p ) is the 
plateau potassium current and J& = Gb(V — Et>) is the 
background current. In these expression, m, h, j, d, f and 
x are gating variables describing the opening and closing 
of ionic channels and the dynamics of these variables is 
described by nonlinear ordinary differential equations of 
the form: 



dy 
dt 



yoo(V) - y 

Ty{V) 



(A3) 



where y represents one of the gating variable. Finally, 
the equations for the currents are supplemented by an 
expression for the calcium concentration: 



[Ca]i 
dt 



-10 4 J si +0.07(10- 4 [Ca] l ) 



(A4) 



Explicit expressions for the constants can be found in 
the original work of Luo and Rudy (Ref. \^)- Here, 
we have made two modifications to allow propagation of 
spiral waves in a small system. First, we have reduced 
G S i from 0.09 to 0.055. Second, we have sped up the 
calcium dynamics by altering the time scales of the d 
and / gates. Specifically, both and Tf were multiplied 
by 0.8. 

Initial conditions that were used consisted of a cable 
initially polarized that was then paced with decreasing 
intervals between stimuli until the desired period was 
reached. 
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